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Abstract 


A  first  order  diffraction  analysis  of  an  optical  interferometer,  Coherent  Gradient  Sensor 
(CGS),  for  measuring  surface  gradients  is  presented.  Its  applicability  in  the  field  of 
fracture  mechanics  is  demonstrated  by  quantitatively  measuring  the  gradients  of  out-of- 
plane  displacements  aroimd  a  crack  tip  in  a  three  point  bend  fracture  specimen  under 
static  loading.  This  method  has  potential  for  the  study  of  deformation  fields  near  a 
quasi-statically  or  dynamically  growing  crack. 
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1.  Introduction 


Presently,  several  optical  methods  are  being  used  for  measuring  elastic  or  plastic 
crack  tip  deformations  under  quasi-static  or  dynamic  loading  conditions.  Among  the 
commonly  used  techniques,  photoelasticity  [I]  measures  princip2d  stress  differences  while 
moire  methods  and  interferometry  [^-4]  map  in-plane  or  out-of-plane  displacements.  In 
the  method  of  caustics,  non-imiform  displacement  gradients  due  to  crack  tip  deforma¬ 
tions  result  in  the  formation  of  a  shadow  spot  [5,6].  Subsequent  interpretation  of  these 
optical  measurements  through  an  assumed  asymptotic  field  description  enables  one  to 
evaluate  the  so-called  stress  intensity  factor  (SIF)  Ki  or  the  J-integral  which  are  widely 
used  in  characterizing  fracture  behavior  of  materials. 

V 

In  most  of  the  above  methods  [1-4],  one  often  encounters  the  difficult  question  of 
whether  the  chosen  experimental  technique  provides  adequate  control  over  the  sensitiv¬ 
ity  of  measurement.  This  becomes  relevant  because  of  the  wide  range  of  magnitudes 
of  deformation  that  may  occur  near  a  crack  tip.  Typically,  interferometric  methods  are 
preferred  for  the  measurement  of  elastic  deformations  while  for  larger  deformations  the 
resulting  fringe  density  often  overwhelms  the  recording  capabilities.  As  a  result,  geomet¬ 
ric  moire  methods  are  used  when  large  deformations  exist.  Howev^)^  fracture  studies 
one  often  needs  a  method  which  can  satisfactorily  perform  in  both  regimes. 

In  this  paper,  we  propose  an  optical  interferometer  -  coherent  gradient  sensor  (CGS) 
to  measure  in-plane  gradients  of  out-of-plane  surface  displacements  around  a  crack  tip. 
The  method  produces  high  contrast  fringes  aind  provides  some  degree  of  control  on  the 
sensitivity  of  measurement  during  quasi-static  experiments.  In  addition,  it  involves  a 
simple  optical  set-up  and,  when  compared  to  other  interferometric  techniques,  it  is  rel¬ 
atively  insensitive  to  vibrations  which  makes  it  a  potential  candidate  for  dynamic  crack 
growth  applications.  Finally,  the  insensitivity  of  this  method  to  rigid  body  motions  is 
highly  attractive  for  solid  mechamics  applications. 

CGS  employs  the  basic  principles  of  the  so-called  ‘moire  deflectometry’  used  widely 


in  a  variety  of  problems  by  Kafiri  and  his  associates  [7,8].  However,  CGS  takes  advantage 
of  coherent  optics  by  using  wave  front  division  (by  means  of  high  density  Ronchi  rulings) 
and  their  subsequent  interference.  By  incorporating  an  online  spatial  filtering  procedure, 
one  could  not  only  realize  high  contrast  fringes  in  real  time,  but  also  use  high  density 
diffraction  gratings  for  the  purpose. 

In  the  following,  we  present  a  diffraction  analysis  of  the  method  and  experimental 
evidence  to  demonstrate  its  applicability  for  elastic  and  elaato-plastic  crack  tip  deforma¬ 
tion  studies. 

2.  Experimental  Method 
2.1  Experimental  Set-up 

In  Fig.(l)  the  schematic  of  the  experimental  set  up  is  shown.  It  consists  of  a  spec¬ 
ularly  reflective  fracture  specimen  illuminated  by  a  collimated  bundle  of  coherent  laser 
light.  Normal  incidence  is  achieved  using  a  beam  splitter.  The  reflected  object  wave  is 
incident  on  a  pair  of  high  density  Ronchi  rulings,  Gj  and  Gj,  separated  by  a  distance  A. 
The  principal  direction  of  the  rulings  is  coincident  with  either  the  Xi—  or  Z2*coordinate 
axes  to  obtain  {duzidxx)  or  (^ua/dxj),  respectively,  where  U3(xi,X2)  denotes  the  out- 
of-plane  surface  displacements.  As  shown  in  the  schematic,  the  origin  of  the  coordinate 
system  is  located  at  the  crack  tip  of  the  specimen.  The  field  distribution  on  the  G2  plane 
is  spatially  filtered  by  the  filtering  lens  Li  and  its  frequency  content  is  displayed  on  its 
back  focal  plane.  By  locating  a  filtering  aperture  around  either  the  ±1  diffraction  orders, 
the  information  regarding  the  displacement  gradients  is  obtained  on  the  image  plane  of 
the  lens  L-i-  It  should  be  pointed  out  that  by  using  cross  gratings  instead  of  line  grat¬ 
ings,  one  could  obtain  both  gradient  fields  simultaneously  by  filtering  the  corresponding 
diffraction  orders  on  the  filter  plane.  In  the  following  sections,  a  first  order  diffraction 
amalysis  is  performed  to  demonstrate  that  the  information  displayed  on  the  image  plane 
indeed  corresponds  to  gradients  of  the  out-of-plane  displ2w:ement  u^. 
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2.2  Principle 


Figure  (2)  explains  the  working  principle  of  the  CGS  interferometer.  For  the  sake  of 
simplicity,  and  without  losing  generality,  the  line  gratings  are  assumed  to  have  a  sinu¬ 
soidal  transmittance.  For  the  case  of  a  plane  wave  reflected  from  the  flat  surface  of 
an  undeformed  specimen  and  propagating  along  the  optical  axis,  the  incident  wave  is 
diffracted  into  three  plane  wave  fronts  Eq^Ei  and  E^i  by  the  first  grating  Gi.  The 
magnitude  of  the  angle  between  the  propagation  direction  of  Eq  and  E±i  is  given  by  the 
grating  equation  6  =  sin~^(A/p),  where  A  is  the  wave  length  and  p  is  the  grating  pitch. 
Upon  incidence  on  the  second  grating  G2,  the  wave  fronts  are  further  diffracted  into  Eq^q, 
Eo,i,  Ei,-i,  Eifi,  £1,1  etc.  These  wave  fronts  which  are  propagating  in  distinctly  different 
directions,  are  then  brought  to  focus  at  spatially  separated  diffraction  spots  on  the  back 
focal  plane  of  the  filtering  lens.  The  spacing  between  these  diffraction  spots  is  directly 
proportional  to  sin  6  or  inversely  proportional  to  the  grating  pitch  p. 

Now,  consider  a  plane  wave  normally  incident  on  a  deformed  specimen  surface.  The 
reflected  light  bundle  incident  on  G\  now  carries  surface  displacement  gradient  informa¬ 
tion,  and  is  constituted  of  light  rays  travelling  in  arbitrary  directions.  If  a  lairge  portion 
cf  such  a  bundle  of  light  has  its  light  rays  nearly  parallel  to  the  opticzJ  axis,  each  of  the 
diffraction  spots  will  be  locally  surrounded  by  a  dispersed  light  field  due  to  the  deflected 
rays.  The  extent  of  this  dispersion  depends  on  the  angle  of  deflection  of  the  reflected  rays. 
By  using  a  two  dimensionaJ  aperture  at  the  filtering  plane,  information  existing  around 
one  of  the  spots  can  be  further  imaged.  This  leads  to  an  important  but  subtle  point 
that  should  be  noted.  Since  each  of  the  diffraction  spots  is  surrounded  by  dispersed  light 
contauning  surfaw:e  deflection  information,  overlapping  of  the  information  could  occur  on 
the  filtering  plane  when  the  deflection  of  the  ray  is  sufficiently  large  (  i.e.,  >  (A/2p)). 
However,  as  will  be  shown  in  the  following  sections,  this  limitation  can  easily  be  overcome 
bv  the  use  of  higher  density  gratings. 
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2.3  Analysis 


Consider  a  specimen  whose  reflective  surface  occupies  the  (xijij)  plane  in  the  unde¬ 
formed  state.  Upon  deformation  the  reflector  surface  C2m  be  «q}ressed  as, 

F{Xi,X2,X3)=X3  +  f{xi,X2)  =  0.  (1) 

Consider  now,  a  plane  wave  which  is  incident  on  the  specimen  along  the  —X3  direction. 
The  unit  surface  normal  N  at  a  generic  point  0(xi,X2)  is  given  by, 

where  e,-  denote  the  unit  vectors  of  the  Cartesian  coordinate  system  (see  Fig. (3))  and 
implies  differentiation  with  respect  to  x^.  Let  d  be  the  unit  vector  along  the  reflected 
ray  whose  direction  cosines  are  cco,  0o  and  70.  From  the  law  of  reflection,  the  coplanar 
unit  vectors  d,  N  and  63  are  related  by  d  •  N  =  63  •  N.  This  leads  to 


d  =  (2e3-N)N-e3. 


(3) 


By  substituting  Eqn.(2)  in  (3), 


d  =  (aoCi  -I-  ^ej  +  7063)  = 


2(/.iei  -h  /,2e2  +  63) 


-63. 


(1  +  /.I  +  /.!) 

Thus,  the  direction  cosines  of  d  can  be  related  to  the  gradients  of  the  function  /  by 

2/.1 _  ,  _  2/3 


(4) 


Oo  = 


(i+fi+fir 


00  = 


(i  +  /.i  +  /S)’  ^  (1  +  /3  +  /J)' 


(5) 


The  ray  along  d  upon  incidence  on  the  grating  Gi  is  split  into  rays  propagating 
along  do,  d±i  whose  amplitudes  £^o(x'),  Ei{x')  and  E-i{xf)  can  be  represented  by, 

iSo(x')  =  aoexp[tido  •  x'],  Ei{x')  =  exp[tfcdi  •  x^,  i?-i(x')  =  oi  exp[tfcd_i  •  x'],  (6) 
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where  cq  and  oj  are  constants  and  the  wave  number  k  =  27r/A.  Due  to  diffraction  by  the 
sinusoidal  grating  Gi,  the  propagation  directions  of  the  resulting  wave  fronts  are  related 
by, 

d±i  =  n^ido,  (7) 

where  are  rotation  tensors  whose  matrices  of  components  axe  given  by, 


=  0  cos^  ±sin^  ,  ( 

^0  ^sin$  cos  6  ^ 

and  6  =  sin~^(A/p).  From  Eqns.(7)  and  (4)  we  find 

cl±i  =  [ofoei  +  cos  5  ±  7o  sin  d)e2  +  (70  cos  0  sin  5)03] .  ( 

Now,  with  reference  to  Fig.(3),  on  plane  G2  we  have 

Eq(O'B)  =  aoexp[tfcdo  •  O'B],  (1 

In  addition,  since  \0'B\do  •  63  =  \0'B\‘yQ  =  A,  one  finds  that  \0'B\  =  (A/70).  Hence, 


Eq{0'B)  =  ooexp  tT; 


’*K). 


In  like  manner. 


|0'i4|di  •  63  =  10'A|(7o cos9  —  ^ sin 0)  =  A, 
10'fi'|d_i  •  63  =  10'S'|(7o  cos  ^  +  /3o  sin  9)  =  A, 


and  thus. 


E\{0'A)  =  Cl  exp[t7;di  •  O' A]  =  oj  exp 


(7ocos^  —  ^sin^)J  ’ 


E^i{0'B')  =  ai  exp[iArd_i  •  0'B'\  =  ai  exp  ik 


(70  cos9  +  j3o  sin  9) 
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2.4  Spatial  Filtering 


Let  the  field  distribution  on  the  G2  plane  due  to  the  amplitudes  Eq,  E\  and  E-i 
be  where  the  double  primed  coordinate  system  is  obtained  by  the  translation 

of  (xi,Z2)  coordinate  system.  Further,  let  the  transmission  function  of  the  sinusoidal 
grating  Gj,  whose  principal  direction  is  parallel  to  the  Xj-axis,  be 

,  27rx^\ 

1  +  cos— (16) 

where  6  is  a  constant.  Then  the  field  distribution  immediately  after  G2  is 

R{x'(,x'')  =  6  +  cos  ^  j  [Eo  +  El  +  E.i]  =  b(l+  cos  ^(xi',x").  (17) 

The  lens  Lx  performs  a  Fourier  transform  of  the  above  field  distribution.  Thus,  the 
amplitude  distribution  on  the  filtering  plane  is 


=  b  -  i)  +  +  i)]  (18) 

2  P  2  p  J 

where  w,-  =  x"/ Xfu  denote  the  in-plane  coordinates  of  the  filtering  plane,  fn  is  the  focal 
length  of  the  filtering  lens  Lj,  5  is  the  Dirac  function,  and  ★  denotes  convolution.  It  is 
clear  from  the  above  expression  that  the  information  regarding  the  surface  deflections 
(direction  cosines  of  d),  which  is  contained  in  H{xx,  Xj),  can  be  filtered  out  by  locating 
suitably  an  aperture  centered  at  =  ±(l/p)  on  the  filtering  plane.  Imaging  this  infor¬ 
mation  through  the  second  lens  £3  is  equivalent  to  performing  a  Fourier  transform  on 
the  filtered  field  which  gives  us  the  amplitude  distribution  on  the  image  plane  as, 

=1  H{-x';,-x';)exp^±i?^y  (w) 

The  intensity  distribution  on  the  image  plsuie  can  now  be  calculated  as, 


/im  =  R^imEin 


H  +  2ii  +  26.6.  CO,  {*A 

I  L  7o(7o  cos  ^  -  A)  sin  0) 

+26,6,  CO,  ( kA  1)  I  ] 

I  [  7o(7oCOS^ -f  ^sin^)  J  J 

+26;  CO,  ( kA  ,J  I . 

I  L  (7o  ^-00  0)  f 
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where  6o  and  bi  axe  constants  and  Ei^  is  the  complex  conjugate  of  Eim-  Under  small  6 
approximation  the  above  equation  reduces  to 


/-  =  ii  +  2i;  +  CO,  cos  • 


(21) 


Thus,  lim  denotes  an  intensity  variation  on  the  image  plane  whose  msodma  occur  when 

kA0o0 


76 


=  2nir, 


n  =  0,  ±1,  ±2, .... 


(22) 


From  Eqn.(5)  the  direction  cosines  qq,  ^  and  70  axe  related  to  the  gradients  of  the 
function  /.  Then,  the  above  equation  can  be  rewritten  as, 

+  |v/P 


kA9 


iv/pp 


=  2nT. 


(23) 


When  [V/p  <<  1,  Eqn.(23)  becomes. 


'-(S 


n  =  0,  ±1,  ±2, .... 


(24) 


where  the  fact  that  $  «  (A/p)  and  k  =  25r/A  are  made  use  of.  Similarly,  when  the 
gratings  are  oriented  with  their  principal  direction  coinciding  with  the  xi-axis, 

/,l»(^)  m  =  0,±l,±2 .  (25) 

Thus,  Eqns.(24)  and  (25)  are  the  governing  equations  for  the  method  of  CGS.  It  should 
be  noted  that  these  equations  are  similar  in  form  to  those  of  moire  deflectometry  [S]  and 
reflection  moire  methods  [9]  based  on  geometric  optics.  It  is  cleair  from  the  above  two 
equations  that  the  sensitivity  of  the  method  could  be  increased  by  either  increasing  the 
grating  sepairation  distance  A  or  decreasing  the  grating  pitch  p.  However,  for  small  A, 
one  needs  to  exercise  caution  in  order  to  avoid  the  formation  of  gap  moire  fringes  when 
the  lens  system  is  arramged  such  that  the  object  surface  is  in  focus.  Another  interest¬ 
ing  point  to  be  noted  is  that  Eqn.(21)  provides  an  explanation  for  the  high  contrast  of 
the  fringes  obtained  in  this  method.  Under  the  assumption  that  the  gratings  used  have 
sinusoidal  transmittance,  the  expression  for  intensity  distribution  /,>„  indicates  a  fringe 
sharpening  effect  due  to  the  existence  of  higher  frequency  term  (2A^^/7o).  The  use  of 
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Ronchi  gratings  with  nearly  rectangular  transmittance  leads  to  further  fringe  sharpening 
due  to  the  existence  of  several  of  such  higher  order  frequencies. 

3.  Experimental  Results 


Two  experiments  are  performed  to  demonstrate  the  applicability  of  CGS  in  solid 
mechanics  in  general  and  fracture  mechanics  in  particular.  First,  the  gradients  of  a 
known  function  /  were  used  to  test  the  reliability  of  the  measurements.  For  this  purpose, 
a  spherical  wave  front  was  generated  using  a  convex  lens  of  focal  length  //=546  mm.  Line 
gratings  of  density  40  lines  per  nun  were  oriented  with  their  principal  direction  coincid¬ 
ing  with  the  i2-zocis  and  were  separated  by  a  distance  of  A=22inm  to  produce  fringes 
that  are  shown  in  Fig.(4).  The  fringe  spacing  {xijn)  corresponding  to  the  gradients  in 
the  i2—  direction  namely,  /,2  measured  from  the  fringe  pattern  is  0.66  mm/fringe.  The 
spherical  wave  front  emerging  from  the  convex  lens  can  be  described  by, 


and  hence  (in  transmission). 


f{XuX2) 


(26) 


,3.) 

which  corresponds  to  straight  line  fringes  as  shown  in  the  figure.  For  the  experimental 
parameters  used  in  this  demonstration,  the  fringe  density  expected  from  the  above  equa¬ 
tions  is  0.63  mm/fringe  which  is  in  good  agreement  with  the  experimental  measurements. 

Secondly,  CGS  was  used  to  obtain  gradients  of  the  o  *t-of-plane  displacement  uz{xi ,  X2) 
(=  X2))  around  a  deformed  crack  tip  in  a  three  point  bend  fracture  specimen  made 

of  AISI 4340  steel.  The  dimensions  of  the  specimen  and  the  experimental  set  up  are  shown 
in  Figs.(5)  and  (6).  The  specimen  has  an  electro-discharge  machined  notch  which  is  25 
fim  wide  and  30  mm  deep  through  a  10  mm  thick  plate.  The  specimen  was  heat  treated 
to  have  a  yield  stress  of  1350  MPa.  One  of  the  surfaces  of  the  specimen  was  lapped  and 
polished  to  obtain  a  flat  reflective  test  surface.  The  specimen  was  statically  loaded  in  a 
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hydraulic  loading  fraune  in  displacement  control  mode.  The  load  and  the  load  point  dis¬ 
placements  were  measured  during  the  test.  The  optical  set  up  has  two  Ronchi  rulings  of 
40  lines  per  mm  density  amd  aire  separated  by  a  distance  (A)  of  21  mm  giving  a  sensitivity 
of  measurement  6.05x10"'*  rad/fringe.  In  Fig.(7)  firinge  patterns  representing  contours 
of  constant  {duzldx\)  and  (dua/Szj)  corresponding  to  three  load  levels  PIPo=0.38,  0.61, 
0.71  axe  shown  where  P<,  is  the  plane  stress  limit  load  for  the  specimen.  For  the  lowest 
load  level  of  the  three  patterns  shown  (Fig.(7a)),  the  neax  tip  out-of-plame  displacement 
field  was  assumed  to  be  given  by  the  linear  elastic,  asymptotic,  plane  stress  expression, 

vh 

f{r,  =  -U3(r,  ^)  =  ■^  kn(r,  (f>)  +  (r22{r,  (f>)] 

«  cos(<?i/2)  -I-  -t-  o(l),  as  r  ^  0,  (28) 

where  r  =  \Jx\-\-  x\,  <f>  =  tan"*(x2/ii),  h  is  the  specimen  thickness,  E  is  Young’s 
modulus,  V  is  Poisson’s  ratio,  is  the  constant  term  in  the  asymptotic  expansion  for 
the  stresses,  and  Ki  is  the  mode-I  stress  intensity  factor  to  be  determined  by  CGS. 


By  using  Eqns.(24)  and  (25),  the  stress  intensity  factor  can  now  be  obtained  from  the 
fringe  pattern,  Fig.(7a),  through; 

iK,)cas  =  ■  (^)  ■  m  =  0,±l,±2,....  (29) 

In  the  above  expression,  m  is  the  fringe  order  and  {Ki)cgs  denotes  the  experimentally 


m  =  0,  il,  ±2, .... 


obtained  value  for  the  stress  intensity  factor.  Note  that  the  constant  term  (7°i  does  not 
appear  in  Eqn.(29)  because  the  measurements  are  sensitive  to  surface  gradients  only. 


{^i)cGS  can  now  be  computed  from  different  fringes  corresponding  to  different  (r,<^) 
pjurs.  If  the  out-of-plane  displacement  field  of  the  specimen  surface  is  indeed  given  by 
Eqn.(28),  the  experimentally  obtained  {Ki)cgs  should  be  independent  of  the  location 
of  measurement  (choice  of  r,4>  pairs  in  Eqn.(29)).  In  such  a  case,  {Ki)cgs  should  agree 
with  the  stress  intensity  factor  (A’/)2D  obtained  from  boundary  value  measurements  by 
means  of  a  two  dimensional  analysis  [10]. 

However,  near  the  crack  tip  three  dimensional  effects  axe  expected  to  violate  the  plane 
stress  assumption  which  leads  to  Eqn.(28).  Indeed,  the  experimental  results  reported  by 


Rosakis  and  Ravi-Chandar  [12]  and  supported  by  the  analytical  investigations  of  Yang 
and  Freund  [11]  have  shown  that  Eqn.(28)  is  a  good  description  of  the  surface  out-of- plane 
displacements  only  for  radial  distances  greater  than  0.5h.  This  behavior  is  also  evident 
from  the  current  experimental  results  displayed  in  Fig.(8)-  Here,  {Ki)2d  is  obtained 
through  the  boimdziry  load  and  the  results  compiled  m  Ref.[10].  (Ki)cgs  Is  obtained 
from  fringes  at  different  radial  distances  along  the  ^  =  0  line  by  using  Eqn.(29).  It  is 
clear  from  the  figure  that  the  measured  value  of  (Ki)cgs  is  in  good  agreement  with  {Ki)2D 
for  (r/h)  >  0.5.  However,  for  (r/h)  <  0.5,  a  marked  under  estimation  of  the  inferred 
stress  intensity  factor  is  observed.  This  is  consistent  with  the  results  shown  in  Figs. (5) 
and  (6)  of  Ref.  [12]  which  are  obtained  by  using  the  optical  method  of  caustics.  This  is 
to  be  expected  since  both  caustics  and  CGS  rely  on  surface  out-of-plane  displacement 
gradients. 

The  fringe  patterns  in  Figs. (7b, c)  correspond  to  substantial  near  tip  plastic  deforma¬ 
tions  and  as  such  the  use  of  Eqns.(28)  and  (29)  is  inappropriate.  Direct  comparison  of 
the  experimentally  obtained  displacement  fields  with  a  three  dimensional  elastic-plastic 
finite  element  analysis  of  the  specimen  is  underway. 

4.  Conclusions 

This  paper  provides  a  detailed  diffraction  analysis  of  the  proposed  interferometer, 
Coherent  Gradient  Sensor.  The  analysis  indicates  that  the  resulting  interference  fringes 
represent  contours  of  constant  surface  gradients.  Although  the  method  bears  similarities 
in  its  operating  principle  with  moire  deflectometry  [8],  it  has  additional  advantages  be¬ 
cause  it  uses  coherent  optics.  Unlike  moire  deflectometry,  which  is  limited  by  diffraction 
effects,  this  method  utilizes  such  effects  by  making  use  of  high  density  gratings  and  spa- 
tiad  filtering  procedure.  Also,  spatial  filtering  enhances  the  contrast  and  sharpness  of  the 
firinges. 

The  feasibility  of  CGS  as  a  means  of  quantifying  crack  tip  deformation  fields  has 
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been  demonstrated  in  the  elastic  and  plastic  regimes.  Simplicity  of  the  optical  set  up, 
sharp  2Lnd  high  contrast  fringes  and  limited  loss  of  light  intensity  make  CGS  a  suitable 
candidate  for  dynamic  crack  initiation  and  propagation  studies.  Its  relative  insensitivity 
to  laboratory  vibrations  comes  as  an  additional  advantage. 
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Figure  Captions 

Figure  1.  Schematic  of  the  Experimental  Set-up  for  CGS 

Figure  2.  Working  Principle  for  CGS 

Figure  3.  Diffraction  of  a  Generic  Ray  in  the  CGS 

Figure  4.  Contours  of  Constant  for  a  Spherical  Wave  Front 

Figure  5.  Specimen  Geometry  of  the  3-point  Bend  Fracture  Specimen 

Figure  6.  Experimental  Arrangement  for  Fracture  Experiment  Using  CGS 

Figure  7a.  Contours  of  {duzjdxx)  and  {du^fdxi)  for  (P/Po)=0.38 

Figure  7b.  Contours  of  {duzldxx)  and  {duzldx^)  for  (P//())=0.61 

Figure  7c.  Contours  of  {duz/dxx)  and  {du^ldx^)  for  {P/Po)=0.71 

Figure  8.  Experimental  Results  from  {duj/dxi)  fringes  along  (r,  ^  =  0)  line  for  (P/Po)=0.38 
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Figure  1.  Schematic  of  the  Experimental  Set-up  for  CGS 


Figure  3.  Diffraction  of  a  Generic  Ray  in  the  CGS 


Figure  4.  Contours  of  Constant  2  for  a  Spherical  Wave  Front 


Thickness,  h=I.O 


Figure  5.  Specimen  Geometry  of  the  3-point  Bend  Fracture  Specimen 
All  dimensions  are  in  cm. 
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Figure  6.  Experimental  Arrangement  for  Fracture  Experiment  Using  CGS 


Figure  7a.  Contours  of  (du^fdx^)  and  {du3fdx2)  for  {P/ Pq)=0.38 
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Figure  8.  Experimental  Results  from  {du^/dxi)  fringes  along  (r, =  0)  line  for  (P/ /o)  = 
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